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Abstract 

It has been recently understood how to deal with high-energy scattering 
beyond the mean field approximation. We review some of the main 
q ! steps of this theoretical progress, like the role of Lorentz invariance and 

O I unitarity requirements, the importance of discreteness and fluctuations 

of gluon numbers (Pomeron loops), the high-energy QCD/statistical 
physics correspondence and the consequences for the saturation scale, 
^ \ the scattering amplitude and other, also measurable, quantities. 

CO 

The high-energy scattering of a dipole off a nucleus/hadron in the mean field approximation is 
decribed by the BK-equation [1]. The main results following from the BK-equation are the so- 
called geometric scaling behaviour of the scattering amplitude [2-4] and the, roughly, powerlike 
energy dependence of the saturation scale [3,4] which are supported by the HERA data [5,6]. 

Over the last few years, we have had real breakthroughs in our understanding of high- 
energy scattering near the unitarity limit. Namely, we have understood how to deal with small- x 
dynamics at high energy beyond the mean field approximation, i.e., beyond the BK-equation. 
In this work, after briefly introducing the known dynamics in the mean-field case, we discuss 
the main steps of the recent theoretical progress as follows: We start with a discussion of the 
first step beyond the mean field approximation, which was done in Ref. [7] by enforcing the 
BFKL evolution in the scattering process to satisfy natural requirements as unitarity limits and 
Lorentz invariance. The consequence was a correction to the saturation scale and the breaking 
of the geometric scaling at high energies. Then, we explain the relation between high-energy 
QCD and statistical physics found in Ref. [8] which has clarified the physical picture of, and 
the way to deal with, the dynamics beyond the BK-equation. We explain that gluon number 
fluctuations from one scattering event to another and the discreteness of gluon numbers, both 
ignored in the BK evolution and also in the Balitsky-JIMWLK equations [9], lead to the breaking 
of the geometric scaling and to the correction to the saturation scale, respectively. In a next step 
we show the new evolution equations, the so-called Pomeron loop equations [10-12], which 
include a new element in the evolution, the Pomeron loop. Finally, we discuss the possibility of 
phenomenological implications [13-17] of the recent theoretical advances. (For further studies 
on the recent theoretical advances (not discussed here) see also [17-27].) 

1.1 Mean field approximation 

Consider the high-energy scattering of a dipole of transverse size r = (x — y) off a target 
(hadron, nucleus) at rapidity Y = ln(l/x). The Y-dependence of the T-matrix in the mean field 
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Fig. 1: The diagrammatic representation of the BK-equation for dipole-hadron scattering. 



approximation is given by the BK-equation (Y -dependence is suppressed for simplicity) 

J d 2 z Af xyz [-T xy + T xz + T zy - T xz T zy ] . (1) 



<9T xy _ a s N c 
dY ~ 2vr 2 



This equation can be interpreted as follows; If increasing the rapidity of the dipole by dY while 
keeping the rapidity of the target fixed, the probability for the dipole to emit a gluon increases. In 
the large- N c limit the initial quark-antiquark state plus the emitted gluon can be viewed as two 
dipoles - one of the dipoles consists of the inital quark and the aniquark part of the gluon whike 
the other dipoles is given by the quark part of the gluon and the inital antiquark. The probability 
for the spliting of the inital dipole (x — y ) into two doughter dipoles with transverse sizes (x — z) 
and (z — y) is given by the weight in Eq.(Q]), a s N c M xyz / (27r 2 )d 2 zdY , where z is the transverse 
size of the emitted gluon and M xyz = (x — y) 2 /[(x — z) 2 (x — z) 2 ] [28]. On the right-hand 
side of Eq. (Q]), the first three terms (first one is virtual) describe the scattering of a single dipole 
with the target whereas the last term gives the simultanous scattering of the two doughter dipoles 
with the target, as shown in Fig.Q] Without the last term, the BK-equation reduces to a linear 
equation, the BFKL equation, which gives the growth of T xy with rapidity, while the nonlinaer 
term, T xz T zy , tames the growth of T xy such that the unitarity limit, T xy < 1, is satisfied. 

One of the main results following from the BK-equation is the geometric scaling behaviour 
of the T-matrix [2-4] in a large kinematical window 

T(r,Y) = T(r 2 Q 2 s (Y)) , (2) 

where Q S (Y) is the so-called saturation momentum defined such that T(r ~ 1/Q S , Y) be of 
0(1). Eq. §2$ means that the T-matrix scales with a single quantity r 2 Q 2 (Y) rather than de- 
pending on r and Y separatelly. This behaviour implies a similar scaling for the DIS cross 
section, a^ p {Y, Q 2 ) = a^* p (Q 2 /Q 2 S {Y)), which is supported by the HERA data [5]. 

Another important result that can be extracted from the BK-equation is the rapidity depen- 
dence of the saturation momentum (leading-^ contribution) [3,4], 



Qt(Y) = Q 2 Exp 



2a s N c x(A ) y 
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(3) 



where x(A) is the BFKL kernel and An. = 0.372. 

The shape of the T-matrix resulting from the BK-equation is preserved in the transition 
region from weak (T ~ a 2 ) to strong (T ~ 1) scattering with increasing Y, showing a "travelling 
wave" behaviour as sketched in Fig|2l on the left hand side. With increasing Y, the saturation 



region at r ^ 1/Q S (Y) where T ~ 1 however widens up, including smaller and smaller dipoles, 
due to the growth of the saturation momentum. As we will see in the next sections, the situation 
changes a lot once gluon number fluctuations are taken into account. 



2 Beyond the mean field approximation 

2.0.1 Lorentz invariance and unitarity requirements 

Let's start with an elementary dipole of size r\ at rapidity y = and evolve it using the BFKL 
evolution up to y = Y. The number density of dipoles of size r 2 at Y in this dipole, n(ri,r 2 , Y), 
obeys a completeness relation 

n( ri ,r 2 ,Y) = J ^n(n,r,Y/2)n{r,r 2 ,Y/2) (4) 

where on the right hand side the rapidity evolution is separated in two successive steps, y = — > 
y = Y/2 -> y = Y. With 

T(ri,r 2 ,Y) ~ co? s r\ n(rt,r 2 ,Y) (5) 
eq.© can be approximately rewritten in terms of the T-matrix as 

j2 r(ri ' ra ' Y) ) - 2^f l dp {h T{ri ^ Y/2) ) (4 T(r ' r2 ' y/2) ) (6) 

where p = ln(rQ/r 2 ). In Ref. [7] it was realized that the above completeness relations, or, 
equivalently, the Lorentz invariance, is satisfied by the BK evolution only by violating unitar- 
ity limits. This can be illustrated as follows: Suppose that r 2 is close to the saturation line, 
r 2 ~ l/Q s (Y), so that the left hand side of Eq.© is large. On the right hand side of Eq.© it 
turns out that T(r\,r, Y/2)/r 2 is typically very small in the region of p which dominates the in- 
tegral. This means that T(r, r 2 , Y/2)/r 2 must be typically very large and must violate unitarity, 
T(r, r 2 , Y/2) » 1, in order to satisfy ©. 

The simple procedure used in Ref. [7] to solve the above problem was to limit the region of 
the p-integration in Eq.© by a boundary p 2 (Y/2) so that T(r, r 2 , Y/2)/r 2 would never violate 
unitarity, or T{r\, r, Y/2)/r 2 would always be larger than a 2 s . The main consequence of this 
procedure, i.e., BK evolution plus boundary correcting it in the weak scattering region, is the 
following scaling behaviour of the T-matrix near the unitarity limit 

and the following energy dependence of the saturation momentum 



Qi(Y) = QZ Exp 

with 
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Fig. 2: Left-hand side: The "travelling wave" behaviour of the solution to the BK-equation. Right-hand side: The 
T-matrix at two different rapidities, Yi and Y%, for different events (thin lines). The thick lines represent the average 
over the events, (T), at the two rapidities, respectivelly. The shape of (T) becomes flatter with rising rapidity. 

Eq.© shows the breaking of the geometric scaling, which was the hallmark of the BK-equation 
shown in Eq.©, and Eq.© shows the correction to the saturation momentum (cf. Eq.©), both 
emerging as a consequence of the evolution beyond the mean field approximation. 

2.0.2 Statistical physics - high density QCD correspondence 

The high energy evolution can be viewed also in another way which is inspired by dynamics 
of reaction-diffusion processes in statistical physics [8]. To show it, let's consider an elemen- 
tary target dipole of size n which evolves from y = up to y = Y and is then probed by 
an elementary dipole of size r, giving the amplitude T[r\, r, Y). It has become clear that the 
evolution of the target dipole is stochastic leading to random dipole number realizations inside 
the target dipole at Y , corresponding to different events in an experiment. The physical ampli- 
tude, T(ri, r, Y), is then given by averaging over all possible dipole number realizations/events, 
T(n,r, Y) = (T(ri,r, Y)), where T{r\,r, Y) is the amplitude for dipole r scattering off a 
particular realization of the evolved target dipole at Y. An illustration is shown in Fig© the 
right-hand side plot, where the T-matrix for different events (thin lines) and the average over all 
events (thick lines), (T), are shown at two different rapidities, respectivelly. 

The mean field description breaks down at low target dipole occupancy due to the discrete- 
ness and the fluctuations of dipole numbers. Because of discreteness the dipole occupancy can 
not be less than one for any dipole size. Taking this fact into account by using the BK equation 
with a cutoff when T becomes of order a 2 s [8], or the occupancy of order one (see Eq.©), leads 
exactly to the same correction for the saturation momentum as given in Eq.©. The latter cutoff 
is essentially the same as, and gives a natural explanation of, the boundary used in Ref. [7] and 
briefly explained in the previous section. 

The dipole number fluctuations in the low dipole occupancy region result in fluctuations 
of the saturation momentum from event to event, with the strength 



extracted from numerical simulations of statistical models. The averaging over all events with 
random saturation momenta, in order to get the physical amplitude, causes the breaking of the 
geometric scaling and replaces it by a new scaling law, the so-called diffusive scaling, in which 



(P 2 s) - (Ps) 2 OC 



a s Y 



(10) 
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case the scattering amplitude is a function of another variable, 



(T(r,Y)) = f 



( 



Hr 2 Q 2 s (Y)) 
^/a s Y/(Ap)* 



) 



(11) 



This equation differs from Eq.O since Eq.O misses dipole number fluctuations. Note that 
because of the geometric scaling violation, the result in Eq.dTTT) changes the shape as the rapidity 
increases, as illustrated in Fig. [2] (right-hand side) by the decreasing slope of the thick line with 
growing rapidity, in contrast to the solution to the BK-equation in Eq.©. 

The statistical physics/high-density QCD correspondence suggests the following picture 
for the wavef unction of a highly evolved hadron which is probed by a dipole of transverse size r: 
As the hadron is boosted to high rapidities the density of gluons inside the hadron grows. Also 
the fluctuation in gluon numbers, which is characterized by the dispersion in Eq. (TTOb . grows 
with rising rapidity. However, as long as a 1 <C 1, which means Y <C Yds — {Ap) 3 /a s , the 
effects of fluctuations can be neglected and the evolution of the hadron is described to a good 
approximation by the BK-equation. Thus, for Y <C Yds, as shown in FigfJl to the left of the 
saturation line, p -C (p s (Y)) = (ln(Q^(Y) Vq)}, is the "saturation region" with the "large-size" 
(small momentum) gluons at a large density, of order l/a s or the T ~ 1, while the shadowed 
region is the transition region from high to low gluon density, or the front of the at T-matrix 
(geometric scaling regime). At higher rapidities, Y S> Yds, where the fluctuations become 
important, the geometric scaling regime is replaced by the diffusive scaling given in Eq. (fTTT t. 
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Fig. 3: The phase diagram of the wavefunction of a highly evolved hadron. 



2.0.3 Pomeron loop equations 

It was always clear that the BK equation does not include fluctuations. However, it took some 
time to realize that also the Balitsky-JIMWLK equations do miss them. It turned out (see first 
Reference in [11]) that the Balitsky-JIMWLK equations do include BFKL evolution, "pomeron 
mergings" but not also "pomeron splittings", which are represented by the three graphs in Fig. [4] 
for two dipoles scattering off a target, respectivelly. After this insight, the so-called Pomeron loop 
equations [10, 11] have been constructed to account for "pomeron splittings" or dipole number 
fluctuations. 

The Pomeron loop equations can be expressed in a Hamiltonian language, in which case 
one extends the JIMWLK-equation [10], or in terms of scattering amplitudes [1 1], in which case 
the Balitsky equations are extended. In order to be close to the BK-equation discussed in sec. II . 1 [ 
we show the Pomeron loops using the scattering amplitude. In the large- N c limit, they can be 
written either as a stochastic equation of Langevin-type, 



<9T xy a s N c 



dY 2tt 2 



( M xyz [ -^Xy "H ^XZ ~\~ T z y T xz T z y\ 



a, a*N r f . |u — v| 



+ ^y^T J A (*> y l u ' z ) (u _ z)2 V V u V v y uv v(u, v, z; Y) (12) 

u,v,z 

or, equivalently, as a hierarchy of coupled equations of averaged amplitudes, where for simplicity 
we show only the first two of them, which read 

d(T xy ) _ a s N c 



dY 2ir 2 



■ J M xyz [-(T xy ) + <T„) + (T zy ) - (T xz T zy )] 



<>UxJ ^'' / - asNc J M xzt [-(T xz T zy ) + (T xt T zy ) + (T tz T zy ) - (T xt T tz T zy )] 
t 

/ M zyt [-(T xz T zy ) + (T xz T zt ) + ( ) — (T xz T zt T ty }] 



dY 2vr 2 

t 

« s iV c 



2vr 2 



'as\ a a N c . 
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J Tl(x, z, z, y|u, v) (T uv ) (13) 



where the noise is non-diagonal (non-Gaussian) in the first two arguments 

(z/(ui,vi,zi;y)i^(u 2 , V2,z 2 ;Y')) = 6 UlV2 5 U2Vl 5 ZlZ2 5 Y y> (14) 
the triple Pomeron kernel [29] reads 

^(xi,yi,x 2 ,y 2 |u, v) = J V^V 2 , [M uvz A(*i, yi|u, z) «4(x 2 , y 2 |z, v)] (15) 

z 

and a 2 A is the amplitude for dipole-dipole scattering in the two-gluon exchange approximation 
and for large-iV c , with 



(a) (b) (c) 

Fig. 4: Two dipoles scattering off a hadron: (a) BFKL interaction, (b) "pomeron mergings", (c) "pomeron splittings". 

In above equations the integrations are always over transverse sizes, J x = d 2 x d 2 y. 

Last term in Eq.([T2l), containing the non-Gaussian noise v, is new as compared with the 
BK-equation and accounts for fluctuations in the dipole numbers or the stochastic nature of the 
evolution in small-x physics. The hierarchy in Eq.([T3l reduces to the BK-equation only in the 
mean field approximation, i.e., when (TT) = (T)(T). The hierarchy in Eq.(fl"3l). as compared 
with the Balitsky-JIMWLK hierarchy, involves in addition to linear BFKL evolution (FiglUa)) 
and pomeron mergings (FiglUb)), also pomeron splittings (FiglUc)), and therefore, in the course 
of the evolution, also pomeron loops. The three pieces of evolution are represented by the linear 
terms, nonlinear terms and the last term on the right-hand side of the second equation in Eq.([T3l. 
respectivelly, which describes the scattering of two dipoles off a target. 

2.0.4 Phenomenology 

It isn't yet clear at which energy fluctuation/Pomeron loop effects start becoming important. The 
results shown in the previous sections, Eq.® and Eq.(Qj]), are valid at asymptotic energies. A 
solution to the evolution equations, which is not yet available because of their complexity, would 
have helped to better understand the subasymptotics. 

Using the statistical physics/high density QCD correspondence, phenomenological conse- 
quences of fluctuations in the fixed coupling case have been studied, for example for DIS and 
diffractive cross sections [15], forward gluon production in hadron-hadron collisions [16] and 
for the nuclear modification factor R p A [13], in case fluctuations become important in the range 
of LHC energies. Recently, in the fixed coupling case, it has been shown that dipole-proton 
scattering amplitudes which include fluctuation effects seem to describe better the HERA data. 
Also the parameters turn out reasonable: The diffusion coefficient D ~ 0.35 (a 2 = D Y) is in 
agreement with numerical simulations of approximations to Pomeron loop equations [19,24] and 
the saturation exponent A ~ 0.2 (Q 2 = (xo/x) x ) is decreased as expected theoretically. On the 
other hand, allowing the coupling to run, however, within a toy model [24] which is supposed to 
mimic the QCD evolution equations with Pomeron loops, it has been argued that gluon number 
fluctuations/pomeron loops can be neglected in the range of HERA and LHC energies. See also 
Refs. [7, 30] for more studies on running coupling plus gluon number fluctuations. 
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